A simple Lattice Model for hysteresis loops with exchange bias. 
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A simple lattice model that allows hysteresis loops with exchange bias to be reproduced is presented. The 
model is based on the metastable Random Field Ising model, driven by an external field, with synchronous local 
relaxation dynamics. The key ingredient of the model is that a certain fraction / of the exchange constants 
between neighbouring spins is enhanced to a very large value Je- The model allows the dependence of several 
properties of the hysteresis loops to be analyzed as a function of different parameters and we have carried out an 
analysis of the first-order reversal curves. 
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The use of zero-temperature lattice models to 
study hysteresis has been very fruitful for the un- 
derstanding of the role played by disorder in de- 
termining loop shape and Barkhaussen noise dis- 
tribution 1 ; . Recently the well known Random 
Field Ising Model (RFIM), has been extended in 
order to account for the existence of exchange 
bias (EB). This property consists in a shift along 
the field axis of the hysteresis loop. EB has 
been experimentally found in many FM/AFM 
bilayer systems 0- The key ingredient of the 
so called Exchange Enhanced RFIM (EE- RFIM) 
[Hj is the existence of a random fraction of the 
nearest-neighbour exchange interactions (of the 
ferromagnetic layer) which is enhanced to a very 
large value. The Hamiltonian of the model on a 
]V = I x L square lattice (that only represents 
the FM layer) can be written as: 
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where Si — ±1 are the spin variables on each 
site, H is the external driving field, hi are 
quenched Gaussian random fields with zero mean 
and variance a 2 . The exchange constants are 
Jij — 1 (energy unit), except for a random frac- 
tion / which has an enhanced value Jij = Je- 
The rate-independent hysteresis properties of this 
model can be simulated by using synchronous me- 
tastable dynamics [Q. The inset in Fig. ^ shows 




Figure 1. A set of FORCS for the EE- RFIM. 
They correspond to averages over 1000 realiza- 
tions with a = 1.0 and /=0.05. The inset shows 
an example of a hysteresis loop with EB obtained 
with the same parameters, reversing the field at 
H= -4 

an example of a hysteresis loop exhibiting EB. 
The reason behind this EB property is that such 
loops do, in fact, correspond to partial loops. The 
symmetric loop is obtained after applying very 
negative (oc Je) external fields. 

In this short paper we present an analysis of 
the First-Order Reversal Curves (FORC) of the 
EE- RFIM corresponding to the case Je = 20 and 
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/ = 0.05 obtained from the simulation of a L = 50 
lattice. Other interesting properties of this model 
have been published elsewhere 

A set of FORC is shown in Fig. ^ In or- 
der to analyze the structure of this set of FORC 
curves, we compute the so called FORC diagram. 
This is obtained by calculating the mixed second 
derivative of the magentization m at field H b on 
a FORC obtained with reversal point H a 0]: 
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The values of p are usually represented on a di- 
agram which is obtained after rotating the co- 
ordinates according to H c = (Hi, — H a )/2 and 
H u = (H a +Hb)/2. The new co-ordinates allow a 
direct comparison with the Preisach distributions 
5 . Fig. [21 shows the central part of the FORC di- 
agram corresponding to the standard RFIM (a) 
and to the EE- RFIM (b) . The comparison reveals 
that the EB property is reflected by a displace- 
ment of the central peak to the bottom. The ver- 
tical width of this peak increases when the width 
a of the random field distribution also becomes 
larger. The shift to the bottom is also greater for 
increasing /. Furthermore, other interesting fea- 
tures appear on the FORC diagram. These can 
be observed on the wider map in Fig. The two 
main characteristics are the occurrence of an os- 
cillation in the H u < region and a second peak 
in the region of large values of H c (~ 17). An 
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Figure 3. FORC diagram corresponding to the 
set of curves in Fig. 1. 

interesting point to make is that the oscillation 
includes a region of negative values of p. Neither 
this negative region nor the asymmetry of Fig. [3] 
can be interpreted within the framework of the 
Preisach model for which p is assumed to be a 
probability density which should always be posi- 
tive and symmetric about the H u = axis. The 
negatives values of p are associated with the fact 
that the slope of the m(H) curves at — 1 < H < 
increases when H a becomes larger. It will be very 
interesting to analyze whether such singular fea- 
tures occur for experimental systems with EB. 
This will reveal the non-Preisach nature of such 
bilayered systems. We acknowledge fruitful com- 
ments from J.Nogues and H.G.Katzgraber. 
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Figure 2. Averaged FORC diagram of the RFIM 
(a) and the EE- RFIM (b) corresponding to a = 2. 
The scale on the right indicates the values of p 
corresponding to different shades of grey. 
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